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ABSTRACT 

PSR B 1259-63 is in a highly eccentric 3.4 year orbit with a Be star and crosses the Be star 
disc twice per orbit, just prior to and just after periastron. Unpulsed radio, X-ray and gamma- 
ray emission observed from the binary system is thought to be due to the collision of pulsar 
wind with the wind of Be star. We present here the results of new XMM-Newton observations 
of the PSR B 1259-63 system during the beginning of 2004 as the pulsar approached the disc 
of Be star. We combine these results with earlier unpublished X-ray data from BeppoSAX and 
' XMM-Newton as well as with ASCA data. The detailed X-ray lightcurve of the system shows 

, that the pulsar passes (twice per orbit) through a well-defined gaussian-profile disk with the 

half-opening angle (projected on the pulsar orbit plane) AO^ 18.5°. The intersection of 
the disk middle plane with the pulsar orbital plane is inclined at ddisk — 70° to the major axis 
of the pulsar orbit. Comparing the X-ray lightcurve to the TeV lightcurve of the the system 
we find that the increase of the TeV flux some 10-100 days after the periastron passage is 
| unambiguously related to the disk passage. At the moment of entrance to the disk the X- 

r/3. • ray photon index hardens from F =* 1 .8 up to T =a 1 .2 before returning to the steeper value 

T > 1 .5. Such behaviour is not easily accounted for by the model in which the X-ray emission 
is synchrotron emission from the shocked pulsar wind. We argue that the observed hardening 
of the X-ray spectrum is due to the inverse Compton or bremsstrahlung emission from 10-100 
/\ ' MeV electrons responsible for the radio synchrotron emission. 
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1 INTRODUCTION 

PSR B 1259-63 is a ~48 ms radio pulsar in a h ighly eccentric 
(e~0.87), 3.4 year orbit with a Be star SS 2883 jjohnston etafl 



1992) . The system 's distance from Earth is d ~ 2 kpc 
Tjohnston et alll994l) . note, however, that d is uncertain by a factor 
of ~2 based on radio measurements (see, e.g., Manchester, Taylor, 
& Lyne 1993). 

Be stars are well-known to be sources of strong highly 
anisotropic matter outflow. Both a dilute polar wind and a denser 
equatorial disc have been invoked to re concile models for infra- 
red, ultra-violet and optical observations JWaters et alll988l) . Tim- 
ing analysis of the PSR B 1259-63 system shows th at the disc 
of Be star is tilted with respect to the orbital plane fWex et all 
Il998t IWangetalJ 120041) . The properties of the radio 
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from the system are very different close to and fa r from the pe- 
riastr on. Rad io observations of t he 1994, 1997 1 Johnstoii^yil 
I l999|) . 2000 jConnors et ai]|2002t) . and the 2004 Ijohnston et a? 
2005) periastron passages show that when the pulsar is far from 
periastron the observed radio emission is comprised entirely in 
highly linearly polarized pulsed emission from the pulsar itself with 
an intensity practically indepe ndent of the pulsar orbital position 
iMcClure-Griffiths et alj[l998T) . But, starting about 100 days be- 
fore periastron, depolarization of the pulsed emission occurs, the 
dispersion measure and the absolute value of the rotation measure 
increase while the flux density decreases. The pulsed radio emis- 
sion then disappears entirely as the pulsar enters the disc and is 
hidden behind it (relative to the observer) approximately 20 days 
before the periastron passage. Shortly before the disc crossing un- 
pulsed radio emission appears and within several days sharply rises 
to a peak which is more than ten times higher than the intensity of 
the pulsed emission far from the periastron. Afterward the unpulsed 
flux slightly decreases, as the pulsar passes through periastron be- 
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fore reaching a second peak just after the pulsar crosses the disc 
for the second time. The unpulsed radio emi ssion is detected until 
at least 100 days after the p eriastron passage. 1 Johnston et all 19991 
12005. IConnors et al.l20o3) . 

The first detection of the X-rays from PSR B 12 59-63 system 
was done by ROSAT in 1992 iCominskv et al.ll994l) . The most re- 
cent published observations of the system in X -rays were carried 
out with the ASCA sa tellite in 1994 and 1995 <Kasr)i et alJll995L 
lHiravamaetai]|l999h . These observations show that X-ray emis- 
sion is approximately twice as high at the time of the disc crossing 
than at periastron. In soft gamma-rays the system was observed 
only after the periastron passage with CGR O/OSSE and INTE- 
GRAL/lSGKl iGrovdl 19951 : IShaw et alj|2004 . The combined soft 
gamma- and X-ray spectrum is consistent with a power law of pho- 
ton index ~ 1.7. This index also coincides with the radio spectral 
index of the unpulsed emission. No pulsed X-ray emission was de- 
tected from the system. With the standard epoch-folding method 
the upper limit to a pulsed component in 2 - 10 Kev energy range 
was estimated to be about 10% of the unpulsed compo nent close to 
the periastron, and to about 40% during the apastron iKaspi et"ai1 
ll995Hffiravama et alll999h. Not e that s earch for the pulsations in 
P - P space allowed iKaspi et alj 11995T) to reduce an upper limit 
close to the periastron to less than 2% under a certain assumptions 
on the shape of pulse profile. 

The 200 4 periastron passage w as observed in TeV gamma- 
rays by HESS lAharonian et alj2005h . TeV lightcurve of the source 
shows an unexpected behaviour. Contrary to the pr edict ion that the 
TeV flux should be maximum at the periastron iKirk et al]|2002h . 
the flux apparently has a local minimum around this moment. 

Several theoretical models were put forward to explain the 
be haviour of the X-ray lightcurve of the system. In the model 
of iTavani & Aronsl ll997h . the enhancement of synchrotron X-ray 
emission during the disk passage is attributed to enhancement of 
magnetic field at the position o f the contact surface of pul s ar and 
stella r winds. In the model of IChernvakova & Illarionovl 1 19991 
l2000l) the enhancement of the inverse Compton X-ray emission dur- 
ing the disk passage is supposed to be due to the effect of the macro- 
scopic mixing of the stellar and pulsar winds in the disk which 
effectively enhances the escape time of high-energy electrons re- 
sponsible for the observed emission. 

Typical energies of electrons respon s ible fo r the X-ray emis- 
sion are E e ~ TeV in [ T avani & Aronsl 1 19971) model while in 
IChernvakova & Illarionovl Jl999h model E. ~ 10 MeV One can, 
in principle, distinguish between the two models by comparing the 
details of radio, X-ray and TeV lightcurve of the system. We have 
organized a cycle of five XMM-Newton observations as the pulsar 
approached the 2004 periastron. In this paper we present the results 
of our observations along with other XMM-Newton observations 
of the system not previously published. For the completeness we 
also include the analysis of the previously unpublished BeppoSAX 
data . The 2004 periastron passage was extensively observe d in ra- 
dio Jjohnston etai]|2005h and TeV lAharonian et alj|2005h bands. 
Combination of the radio, X-ray and TeV data opens up a possibil- 
ity, to study the evolution of the spectrum of PSR B 1259-63 system 
over some 12 decades in energy. 

Our analysis of large X-ray data set enables, for the first time, 
to constrain the geometrical parameters of the disk of Be star, such 
as the opening angle and orientation. Comparing the TeV and X- 
ray lightcurves we find that the increase of the TeV flux some 10 
to 100 days after the periastron passage is unambiguously related 
to the pulsar passage through the Be star disk for the second time. 
This challenges the conventional interpretation of the observed TeV 
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Figure 1. Schematic representation of PSR B1259-63 binary system with 
the locations of the pulsar during XMM-Newton BeppoSAX and ASCA ob- 
servations. Shaded area shows the geometry of the disk inferred from the 
X-ray data. Darker and lighter shading corresponds to 1 and 2 half-opening 
angles of the disk, respectively. Angle 8 is the orbital phase of the pulsar 
referred in the text. 

emission as the inverse Compton emission from the pulsar wind 
electrons. 

Another surprising result of our analysis that is not predicted 
in the theroetical models cited above is the persistently hard X-ray 
spectrum (photon index r m 1.2-1.3) observed during the 10-day 
period when the pulsar enters the disk. 

This paper is organized as follows: in Section 2 we describe 
the details of the XMM-Newton and BeppoSAX data analysis. The 
results are presented in Section 3, and discussed in Section 4. 



2 OBSERVATIONS AND DATA ANALYSIS 

In FigureQwe show a schematic drawing of the pulsar orbit along 
with the approximate location of the pulsar position during XMM- 
Newton, BeppoSAX, and ASCA observations. The X-ray observa- 
tions are more-or-less homogeneously distributed over the whole 
range of pulsar orbital phase < 6 < 360°. The figure also shows 
the location of the Be star disk (parameters of the disk are found 
below from the analysis of the X-ray data). One can see that the set 
of observations denoted as X6 - X10 enables to study the entrance 
to the disk in great details. 

2.1 XMM-Newton observations 

The log of the XMM-Newton data analyzed in this paper is pre- 
sented in Table 1. Observations X6 - X10 were obtained as part 
of the 2004 periastron campaign and we have also used all other 
XMM-Newton observations of the system, which are now public 
(XI - X5). In the Table negative r denotes the number of days 
before the date of 2004 periastron passage (2004 March 7), and 
positive values correspond to the number of days after the 2000 pe- 
riastron (2000 October 17). 9 is the orbital phase counted as shown 
in Fig.Q] 

The XMM-Newton Observation Data Files (ODFs) were ob- 
tained from the on-line Science Archive 1 ; the data were then pro- 
cessed and the event-lists filtered using xmmselect within the Sci- 
ence Analysis Software (sas) v6.0. 1. In all observations the source 
was observed with the MOS1 and MOS2 detectors. The XI - X5 

1 http://xmm.vilspa.esa.es/external/xmm_data_acc/xsa/index.shtml 
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Table 1. Journal of XMM-Newton observations of PSR B 1259-63 



Data 


Date 


MJD 


T 


e 


Exposure 


Set 






(days) 


(deg) 


(ks) 


XI 


2001-01-12 


51921.73 


87.2 


320.8 


11.3 


X2 


2001-07-11 


52101.31 


266.7 


342.9 


11.6 


X3 


2002-07-11 


52467.24 


-604. 1 


0.5 


41.0 


X4 


2003-01-29 


52668.27 


-403.1 


9.31 


11.0 


X5 


2003-07-17 


52837.53 


-233.9 


19.5 


11.0 


X6 


2004-01-24 


53028.79 


-42.6 


57.4 


9.7 


X7 


2004-02-10 


53045.43 


-25.0 


73.4 


5.2 


x«s 


2004-02-16 


53051.39 


-20.0 


83.1 


7.7 


X9 


2004-02-18 


53052.02 


-18.4 


86.5 


5.2 


X10 


2004-02-20 


53055.82 


-15.6 


93.2 


6.9 



observations were done in the Full Frame Mode, while the 2004 
observations were performed in the Small Window Mode, to mini- 
mize the pile-up problems. For the X6 - X10 observations PN data 
are also available. In all observations a medium filter was used. 

The event lists for spectral analysis were extracted from a 45" 
radius circle at the source position for the XI - X5 observations, 
and from a 22.5" radius circle for MOS 1,2 observations in Small 
Window Mode (X6 - X10). For the PN instrument a region of 35" 
around the source position was chosen. Background photons were 
collected from a region located in the vicinity of the source with 
the same size area as the one chosen for the source photons. 

For the spectral analysis, periods of soft proton flares need to 
be filtered out. To exclude them we have built light curves with 
a hundred second binning and excluded all time bins in which the 
count rate above 10 keV was higher than 1 .5 cnt/s. The source coun- 
trate in 1 - 10 keV energy range varied from about 0.07 to 1.2 cts/s 
for MOS instruments, and from 0.4 to 4 cts/s for PN. Data from 
MOS1, MOS2 and, when available, PN detectors were combined 
in the spectral analysis to achieve best statistics. 

2.2 Beppo SAX observations 

Table [5] lists the observations of PSR B 1259-63 made with 
BeppoSAX around the 1997 periastron passage. In the Table t de- 
notes the number of days before or after the 1997 periastron pas- 
sage. 

We used data of MECS instruments aboard BeppoSAX obser- 
vatory which provides the coverage in the energy range ~ 1.5- 11.0 
keV. Reduction of MECS data was done with the help of standard 
tasks of LHEASOFT/FTOOLS 5.2 package according to recom- 
mendations of BeppoSAX Guest Observer Facility 2 . Spectrum of 
the source was obtained extracting events from a circle of 3' radius 
around the source position. Spectrum of the instrument background 
was estimated considering events from a circle of the same radius 
6' away from the source. Typically the source provided several 
hundreds of counts per observation. Response matrices (RMF and 
ARF) of MECS detectors were taken from HEASARC archive 3 . 
Corrections of photons arrival time to the Solar system barycenter 
for subsequent usage of the event lists in timing analysis were done 
with the help of tasks from SAXDAS package 4 . 

We have checked that all the spectra obtained in the S2 - S5 

2 http://heasarc.gsfc.nasa.gov/docs/sax/abc/saxabc/saxabc.html, 

http://www.asdc.asi.it/bepposax/software/cookbook/ 

i http://heasarc.gsfc.nasa.gov/FTP/sax/cal/responses 

4 http://www.asdc.asi.it/bepposax/software/saxdas/index.html 



Table 2. Journal of BeppoSAX observations of PSR B 1259-63 



Data 


Date 


MJD 


T 


6 


Exposure 


Set 






(days) 


(deg) 


(ks) 


SI 


1997-03-22 


50529.63 


-68.31 


44.9 


28.622 


S2 


1997-09-02 


50693.75 


95.81 


322.9 


16.985 


S3 


1997-09-08 


50699.44 


101.50 


324.1 


53.555 


S4 


1997-09-17 


50708.08 


110.14 


325.9 


51.301 


S5 


1997-09-25 


50716.05 


118.11 


327.3 


28.719 
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Figure 2. Contour plot of the XMM-Newton field of view for the X6 ob- 
servation. A total of 10 contours were used with a linear scale between 5 
counts per pixel (outer contour) and 50 counts per pixel (inner contour). 
In this observation 2RXP J130159. 6-635806 was forty times brighter than 
PSR B 1259-63. This figure is taken from Chernyakova et al. (2005). 

observations are consistent with each other, and combined them in 
order to improve statistics. This combined spectrum we denote as 
SC hereafter. 



3 RESULTS 

3.1 Imaging Analysis 

In the XMM-Newton field of view during its PSR B 1259-63 obser- 
vation program, two sources were clearly detected (see Fig. [2] for 
the contour plot of XMM-Newton field of view for the MJD 53029 
(X6) observation). Besides PSR B 1259-63 itself there is second 
source, which we have identified as 2RXP J13 0159.6-635806. De- 
tailed analysis of this source are presented in IChernvakova et alj 
120051) . In this paper it is shown in particular, that 2RXP J130159.6- 
635806 is highly variable, with an intensity during flares which 
is several times larger than the peak intensity of PSR B 1259-63. 
Therefore, although many observations of PSR B 1259-63 have 
been made with RXTE, the fact that it is a non-imaging telescope 
means that it is extremely difficult to determine the flux from 
PSR B 1259-63 alone. We therefore do not use numerous RXTE ob- 
servations of the 2004 periastron passage, presented in lShaw et alj 
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Figure 3. Comparison between the X-ray (top), TeV (middle) and radio (bottom) lightcurves. The 1-10 keV X-ray flux is given in 10~"ergs cirT 2 s~'. 
XMM-Newton observations are marked with triangles, BeppoSAX ones with circles, and ASCA ones with squares. Data for four different periastron passages 
are shown with different colors: red (1994), green (1997), black (2000) and blue (2004). Bottom X axis shows the orbital phase, 8, top X axis shows days from 
periastron, t. 



120041) . Note that despite in the lShaw et ai]<2004 all RXTE obser- 
vations are analysed, the authors were able to use in their analysis 
only data simultaneous to INTEGRAL observations, when it was 
clear that influence of 2RXP J130159.6-635806 is small. 

3.2 The X-ray lightcurve 

The upper panel of Fig. [3] shows the X-ray lightcurve of the 
system. For comparison we also sho w in the same Figure TeV 
lightc urve of 2004 HESS observation iAharonian et alj|2005h and 
radio 1 Johnston et all 1 19991: IConnors et all 120021: Ijohnston et alj 
2005D lightcurves of different years periastron passages. 

The X-ray flux from the source, in agreement with previous 
observations, was observed to be highly variable. The 1-10 keV 
flux is at minimum near apastron (observation X3), and has max- 
imum more than 20 times larger, 18 days before periastron (X9). 
During observations X5, X6, and X7, the intensity of the emission 
increased, as the pulsar approached periastron and the Be star disc. 
Five days after X7 observation the flux increased by a factor of 4. 



During the next two days the intensity continued to increase, and 
reached its peak value at X9. The next observation (X10) was done 
only three days later, and the observed flux was already one and a 
half times lower. 

Comparison of X-ray and radio lightcurves shows that the X- 
ray flux from the system varies in a more regular way than the ra- 
dio flux: the X-ray data points from different periastron passages 
lie more-or-less on the same curve. Rapid grows of the X-ray flux 
found in our XMM-Newton observations X6 - X10 is correlated 
with the rapid growth of the unpulsed radio emission from the sys- 
tem. The growth of radio and X-ray flux at these phases can be 
attributed to the pulsar entering the Be star disk. 

Unfortunately, TeV observations start somewhat later and it is 
not possible to see whether the TeV flux growth during the pre- 
periastron disk crossing. However, simple geometrical argument 
tells that the orbital phase 9 at which the pulsar should enter the 
disk for the second time should be shifted by 180° relative to the 
first entrance. From Fig.[3]one can infer that the first pre-periastron 
entrance falls roughly between the phases 70° < 9 < 110°. Thus, 
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Figure 4. The X-ray (top) and TeV (bottom) flux as a function of the relative 
phase A8 = 8 - ddi sk (see text for the definition of 6m sk ), The curves show a 
fit with a gaussian of the half-width A8jj sk = 18.5°. 



the pulsar has to enter the disk again between the phases 250° < 
8 < 290°. Surprisingly, one can clearly see from the middle panel 
of Fig.|5|that the TeV flux grows in this phase interval. 

The growth rate of the X-ray flux between the phases 70° < 
8 < 1 10° is remarkably close to the decrease rate at the phases 
290° < 8 < 340°. This enables us to make a conjecture that the 
observed decrease of the X-ray flux can be associated to the exit of 
the pulsar from the disk. To test this conjecture we superimpose the 
pre-periastron X-ray and TeV lightcurves over the post-periastron 
lightcurves by shifting the phase of the post-periastron data points 
by -180°. For X-rays we have used only data from the rise (8 < 
110°) and the fall (8 > 290°) periods. The result is shown in Fig. 
One can see that in such representation the rise and decrease of 
both X-ray and TeV flux from the system can be well fitted with 
a gaussian curve F(0) ~ exp {-(8 - 8 disk ) 2 l(2A8 disk )}. We find that 
the best fit is achieved with the parameter choice 8 disk m 109.1°, 
A8 d j sk — 18.5° (the coordinate A8 along the X-axis of Fig.[4]is, in 
fact AO = 8- 8 disk ). 

The most obvious explanation to the shape of the X-ray and 
TeV lightcurves shown in Fig. [4] is that pulsar crosses the thick 
disk with the half opening angle A8 disk . The "middle" plane of 
the disk intersects the pulsar orbital plane along the line 6 = 
(109.1° U 289.1°). The intersection of the disk with the above pa- 
rameters with the pulsar orbital plane is shown schematically in 
Fig. Q by a shaded area. Because of the exponential profile, the 
disk has no fixed boundary. Denser shading in Fig.Qcorresponds 
to the one opening angle of the disk, while lighter shading corre- 
sponds to two opening angles. The disk is also shown by the same 
shading in the Fig. [3] From this figure one can see that although 
the disk appears quite symmetric in terms of the orbital phase 8, it 
is highly asymmetric in terms of the time r which measures days 
from the periastron passage. For example, during the second (post- 
periastron) disk passage, the rise phase (250° < 8 < 290°) of X- 
ray and TeV lightcurves takes just some 10 days, while the decay 
(290° < 8 < 340°) takes about 100 days. 
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Figure 5. Folded (upper panel) and unfolded (bottom panel) PSR B1259- 
63 spectra from the SC, X3, X7, and X9 observations. The XMM-Newton 
spectra for X7 and X9 are for MOS 1 instrument. The X3 spectrum is for 
PN instrument. 



3.3 Timing Analysis 

For each of the five PN event files in 0.2 - 10 keV energy range 
with a 6 ms time resolution we searched for a possible pulsation 
around the nominal 47.6 ms period using the epoch-folding tech- 
nique (efsearch tool from the XR0N0S package). From the indi- 
vidual observations the best 3 - cr upper limit on pulsation at a 
period of 47.6 ms is ~ 9%, compatible with previous similar an aly- 
sis made with ASCA <Kaspi et alll995tlHiravama et alll999h . In 
order to perform a deeper search for pulsations we combined neigh- 
boring XMM observations (namely X7-X10, see. Table 1). As the 
time span of these observations covers only small fraction of the 
binary orbit we did not make corrections to the binary motion of 
the neutron star. Usage of the combined XMM data improves the 
upper limit to ~2%. 

Data of BeppoSAX/MECS (1.5-11 keV) observations also 
were searched for pulsations. For deeper search we combined ob- 
servations S2-S5 (see Table 2). Due to small time period of per- 
formed observations we did not make correction for binary motion 
of the neutron star. In spite of considerably smaller effective area 
of MECS instruments of BeppoSAX observatory, the effective ex- 
posure time of observations of the source was much larger, which 
resulted in a comparable upper limit on X-ray pulsation fraction. 
We have not found pulsations with 3-cr upper limit at ~2%. 



3.4 Spectral Analysis 

The spectral analysis was done with NASA/GSFC XSPEC soft- 
ware package. In Fig. [5] the folded and unfolded spectra of 
PSR B 1259-63 for SC, X3, X6, and X9 observations are shown. 

A simple power law with a photoelectrical absorption de- 
scribes the data well, with no evidence for any line features. In Ta- 
ble|3|we present the results of the three parameter fits to the XMM- 
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Figure 6. Spectral parameters of XMM-Newton, BeppoSAX, and ASCA observations of the PSR B1259-63 system. Notations are the same as in Fig. [3] Top 
panel shows the evolution of the photon index T as a function of the orbital phase 6. Bottom panel shows the evolution of the hydrogen column density Nh- 
Denser and lighter shaded region show phases within one and two half opening angles of the Be star disk. 



Newton data in the 0.5 - 10 keV energy range. The uncertainties 
are given at the lcr statistical level and do not include systematic 
uncertainties. The graphical representation of the evolution of the 
PSR B 1259-63 spectral parameters along the orbit is given in Fig- 
ure [6] which includes data from the BeppoSAX an d XMM-Newton 
observ ations, and those from ASCA (taken from iHiravama et alj 

In Figure^we show contour plots of N H versus photon spec- 
tral index T, holding the 1-10 keV flux fixed at the values given 
in Table[3] The contours plotted are the 68%, 90%, and 99% con- 
fidence levels. For the X5 observation the values of Nh and T are 
consistent with X4 observation, and we omitted the X4 contours to 
keep the Figure clear. 

One can see that the variation of N H as a function of orbital 
phase is correlated with the Be star disk crossing. No variations 
of N H are observed for the XMM-Newton observations X2 to X7 
which all fall outside the disk (assuming the disk parameters found 
in Section 3.2). The mean value of of N H during these observations 
is N H (X2 - XI) ~ 0.32 x 10 22 cnr 2 . However, the observations X8, 
X9, X10 and XI are characterized by absorption which is higher by 
a factor of about 1.5, N H (X% - X10,XV) ~ 0.45 x 10 22 cnT 2 . Thus, 



the hydrogen column density remains high starting from about 20 
days before the periastron till at least 87 days after periastron. 

The behaviour of the spectral slope as a function of the orbital 
phase 6 found in XMM observations differs from the one found 
in ASCA observations. Whereas in 1994 - 1995 the spectrum at 
apastron was found to be significantly harder than during perias- 
tron, XMM-Newton observations show much softer spectrum at the 
apastron with the spectral slope close to the ASCA periastron value. 

The most remarkable feature of the spectral evolution of the 
system is the hardening of the X-ray spectrum close to the moment 
when pulsar enters the Be star disk at the phase 8 = - 2A# AjJ . — 
70° . One can see that the decrease of the photon index T is simul- 
taneous with the onset of the rapid growth of the X-ray flux. To 
the best of our knowledge, such behaviour was not predicted in any 
of existing models of X-ray emission from the system. Moreover, 
the observed values T < 1.5 during the observations X6 - X10 
are difficult to reconcile with the models of synchrotron or inverse 
Compton emission from shock-accelerated electrons because of the 
too hard electron spectrum implied. Possible error in determination 
of T due to the uncertainty in determination of N H (see Fig.0 is 
too small to allow for T > 1.5 at least for observations X7 - X9. 
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Figure 7. Confidence contour plots of the column density Nh and of the photon spectral index I" uncertainties for a power-law fit to XMM-Newton observations. 
The contours give 68%, 90%, and 99% confidence levels. X4 observation, marked with triangle, is very similar to the X5 observation and the contours are not 
plotted here. 



Table 3. Spectral parameters for BeppoSAX and XMM-Newton Observa- 
tions of PSRB 1259-63.* 



Data 


f (1-10 keV) 


r 


N H 


X 1 (dof) 


Set 


l(T 12 erg s _I 




(10 22 era- 2 ) 




SI 


1.36+0.26 


1 oq+0.73 
■ -0.08 


47+1.48 


0.814(201) 


SC 


631 +8:S 

U - J -0.80 


65+" 16 


0.898(813) 


XI 


8.51(30) 


1.51(3) 


0.45(2) 


0.97 (341) 


X2 


2.81(23) 


1.36(7) 


0.39(4) 


1.02(122) 


X3 


1.04(49) 


1.69(4) 


0.29(2) 


0.88 (207) 


X4 


1.82(14) 


1.82(7) 


0.36(3) 


1.01 (109) 


X5 


0.99(9) 


1.80(8) 


0.32(4) 


0.80 ( 62) 


X6 


2.47(17) 


1.41(6) 


0.35(3) 


0.96(151) 


X7 


4.05(23) 


1.20(4) 


0.28(3) 


0.95 (154) 


X8 


12.59(34) 


1.31(2) 


0.47(1) 


1.12(614) 


X9 


24.27(56) 


1.39(2) 


0.43(1) 


1.07 (737) 


X10 


16.79(37) 


1.47(2) 


0.45(1) 


0.88 (772) 



* Number in parentheses represent 68% confidence interval uncertainties in 
the last digit quoted. 



4 DISCUSSION 

In all theoretical models of unpulsed emission from PSR B 1259-63 
the collision of relativistic pulsar wind with non-relativistic wind of 
the Be star plays a crucial role. Due to the interaction of the winds a 
system of two shock waves forms between the stars. Close to the Be 
star and to the pulsar both winds are supposed to be radial but after 
passing the shocks particles turn and start to flow along the contact 
surface, losing their energy via synchrotron and inverse Compton 
(with the seed photo ns being the Be star soft photons) emissio n. 

In the model of IChernvakova & Illarionovl <1999ll200Cl) (CI) 
the observed X-ray emission is due to inverse Compton scattering 
of soft photons from the Be star on the same electron populaiton 
that produce the non-pulsed radio emission from the system. The 
intensity of observed unpulsed X-ray and radio emission depends 
strongly on the flow velocity of the relativistic particles beyond the 
shock. Slower drift velocity leads to the higher concentration of the 
relativistic electrons beyond the shock and higher X-ray and radio 
luminosity. Decrease of the drift velocity during the pulsar passage 



through the disk is caused by the macroscopical mixing of pulsar 
wind with the stellar wind. 

Main prediction of this model is the correlation of X-ray with 
the radio lightcurves of the system. Although different quality of 
X-ray and radio lightcurves does not allow to make definitive con- 
clusions about radio-X-ray correlations, our XMM-Newton obser- 
vations show that the onset of rapid increase of the X-ray flux is 
simultaneous with appearance of the unpulsed radio emission from 
the system (see Fig[3}. 

The X-ray inverse Compton emission in the CI model is pro- 
duced by electrons with rather moderate gamma-factors, y e IC ~ 10 
while the gamma-factors of radio-synchrotron-emitting electrons 
are around y e j ~ 10 2 . The synchrotron/inverse Compton cooling 
times for such electrons are 



tic 



10 s 
10" 



0.1 G 



B 



3/2 



10 GHz 



1 erg/cm 



1 keV 



eic 



1/2 



(1) 



Both are large compared to the typical adiabatic cooling (or escape) 
time t esc ~ R/v ~ 10 3 s (R is the characteristic size of the emission 
region, and v is the velocity of the flow along the contact surface). 
Another important time scale for the electrons in the energy range 
relevant in the CI model is the Coulomb loss scale 



t coulomb ~ 10 



10 9 cm 3 



(2) 



As the pulsar enters the disk, the density of the stellar wind in- 
creases. Typical estimates of the wind denisty at the pulsar location 
(assuming the radial profile n ~ (R+/R) 2 ) result in the density esti- 
mate n ~ 10' - 10 10 cirr 3 . Thus, inside the disk Coulomb cooling 
time is comparable to the escape time. 

This effect leads to the modification of electron spectrum at 
the energies below the so-called "Coulomb break" energy at which 

t coulomb tesc- 



I ejcoulomh 



10" cm 



ir t esc 

>Hl0 3 s 



(3) 



sJ 

Namely, since the Coulomb loss rate dE/dt ~ y e lt coll i omD is energy- 
independent, the electron acceleration spectrum with the spectral 
index Y e > 2 hardens to (T e — 1) < 2 below the break energy. This 
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Figure 8. Modification of the broad-band spectra of CI (a) and TA (b) 
models due to the bremsstrahlung emission from shocked stellar/pulsar 
wind electrons during the disk passage. Contributions to the total spec- 
trum from synchrotron, inverse Compton and bremsstrahlung processes 
are marked with S, IC and BS, respectively. In the case of TA model 
black solid curves show the spectrum of emission from the shocked pul- 
sar wind while the green dashed curves show the contribution from the 
shocked stellar wind. Parameters used for the CI model calculation are 
n = 2 X 10'cnr 3 , B = 0.01 G (low) and B = 0.1 G (high synchrotron 
flux), dN e /dy ~ y~ 2 - 4 exp(— y/400), with such parameters Coulomb break 
is at y ~ 50. The same parameters are used for the calculation of the 
shocked stellar wind emission (green dashed curves) in the lower panel. 
The emission from the shocked pulsar wind was calculated assuming n = 
2 X 10 u cnr 3 , B = 0.15 G (low) and B = 1.5 G (high synchrotron flux), 
dN e ldy ~ y~ 2 - 3 exp(-y/4 X 10 6 ). With such parameters Coulomb break is 
at y ~ 500. 



naturally leads to the hardening of the spectrum of inverse Compton 
X-ray emission at the moment when the pulsar enters the disk, as 
it is shown in Fig. [8^. One can note that the Coulomb break energy 
is close to the energy of the X-ray emitting photons. This means 
that the X-ray photon index is very sensitive to the slight changes of 
the density and/or of the outflow velocity (see Fig.|8^ in which the 
two spectra are calculated for y ef0 uiomb = 10 and y efou i om ), = 50). 
Since both vary during the disk passage (it is assumed that within 
the disk pulsar and stellar winds are macroscopically mixed), nu- 
merical modelling is required to predict the time evolution of the 
X-ray spectr al index. We leave the detailed modelling of this for 
future work IChernvakova & Neronovl200r3) . 

In the original CI model, only electorns with typical gamma- 
factors of 10-100 are considered. The presence of protons with 
comparable gamma-factors y p ~ 1 - 10 3 at the contact surface 
would lead to the generation of TeV emission during the disk pas- 
sage which is due to the proton-proton interactions in the disk (see 
lKawachietai]<2004 ). Macroscopic mixing of the pulsar and stel- 
lar wind will lead to the increase of the proton-proton interaction 
rate in the disk, which can explain the form of the TeV lightcurve 
of the system. 

In the model of iTavani & Aronsl 1 19971) (TA) the X-ray emis- 
sion is the synchrotron emission from shock-accelerated electrons 



of the pulsar wind. Typical gamma-factors of the X-ray emitting 
electrons are y ~ 10 6 . Modulation of the X -ray flux from the source 
during the disk and periastron passage is due to the combined effect 
of displacement of the contact surface by the disk (which leads to 
the variation of the magnetic field strength) and the increase of the 
inverse Compton energy loss close to the periastron. 

One of the main predictions of TA model is steepening of the 
X-ray spectrum during the disk and periastron passage. The effect 
is due to the fact that synchrotron and/or inverse Compton cool- 
ing time of X-ray emitting electron becomes comparable or shorter 
than the time needed to escape from the emission region. However, 
the observed behaviour of the spectrum is quite opposite: the en- 
trance to the disc is accompanied by the hardening of the spectrum. 

The values of the photon index T < 1.5 found in the observa- 
tions X6 - X10 are too hard to be explained with the synchrotron 
emission from shock-accelerated electrons which normally result in 
the electron spectra with power-law index T e > 2 and, correspond- 
ingly in synchrotron spectra with photon index r > 1.5. Of course, 
over a large fraction of the pulsar orbit the synchrotron cooling time 
is larger than the time needed to escape from the emission region 
along the contact surface. Thus, if one assumes that the produc- 
tion spectrum of high-energy electrons responsible for the X-ray 
emission is harder than T e = 2, (if, e.g. electrons originate in the 
cold pulsar wind with the bulk Lorentz factor y ~ 10 s ), it is possi- 
ble, in principle, to obtain spectra which are harder than T = 1.5. 
However, in this case one expects to observe the very hard X-ray 
spectrum with T < 1.5 over most of the pulsar orbit. 

Within the TA model radio synchrotron emission originates 
from different population of electrons than the X-ray synchrotron 
emission (radio is convensionally attributed to shock-accelerated 
stellar wind electrons, while X-ray flux is emitted by the pulsar 
wind electrons). Although it is not possible to explain the harden- 
ing of the X-ray spectral index by the changes in the pulsar wind 
electron spectrum, one can try to ascribe the observed hardening to 
the contribution from the shocked stellar wind electrons. The X-ray 
spectrum with photon index harder than T = 1.5 is natural if there 
is a bremsstrahlung contribution to the X-ray flux. The increase of 
the density of stellar wind inside the disk can lead to the increase 
of bremsstrahlung luminosity. 

The estimate of the bremsstrahlung loss time 



shows that for radio emitting electrons it is orders of magnitude 
shorter than the synchrotron cooling time (see Q). Thus, the 
bremsstrahlung luminosity should be much higher than the radio 
synchrotron luminosity. From Fig.|8|one can see that a satisfactory 
fit to the broad-band spectrum can be obtained assuming the disk 
density n ~ 10" cm" 3 . It is interesting to note that bremsstrahlung 
can also explain the behaviour of the TeV lightcurve in this case, if 
one assumes the possibility of mixing of the pulsar wind electrons 
with the pulsar disk medium (similar to the one assumed in the CI 
model). 

The bremsstrahlung model faces, however, a serious over- 
all energetics problem. From Fig. [SJ> one can see that the total 
bremsstrahlung luminosity is at least 2 orders of magnitude higher 
than the X-ray luminosity at the moment of X7 observation, L X i — 
10 33 erg/s (see Fig. [5). As it is discussed above, the energy loss rate 
of mildly relativistic electrons responsible for the bremsstrahlung 
emission is dominated by the Coulomb losses. Comparing to 
one can see that the Coulom loss rate is still 2-3 orders of mag- 
nitude higher than the bremsstrahlung luminosity. Thus, to obtain 



XMM-Newton observations ofPSR B1259-63 near the 2004 periastron passage 9 



the X-ray bremsstrahlung flux at the level of 10 33 erg/s one has to 
assume that the power injected into the 10 - 100 MeV electrons is 
about 10 38 erg/s which is much higher than the spin-down luminos- 
ity of the pulsar L pll i sar ^ 10 36 erg/s. 



5 CONCLUSIONS 

In this paper we have presented the XMM-Newton observations of 
the pulsar PSR B 1259-63 during its 2004 periastron passage near 
the companion Be star SS 2883. Combining our observations with 
the previous X-ray observations from XMM-Newton, BeppoSAX 
and ASCA we produced the detailed lightcurve of the system over 
the wide range of the orbital phases. 

Using the large X-ray data sample we were able to constrain 
the geometry of the Be star disk. In particular, we have found that 
the half-opening angle of the disk (projected on the pulsar orbital 
plane) is A6,u s k - 18.5° and that the line of intersection of the disk 
with the pulsar orbital plane is inclined at about 70° with respect to 
the major axis of the pulsar orbit. 

Inspection of the X-ray and TeV lightcurves of the system 
has revealed correlation of the X-ray and TeV flux variations. Our 
analysis shows that the variations of the TeV flux some 10-100 
days after the periastron passage can be well fit with the gaussian 
curve whose parameters coincide with the parameters of the disk 
found from the X-ray analysis (center at G disk = 289.1°, half-width 
A9 = 18.5°) (see Fig. |4j- If the observed behaviour of the flux is 
indeed due to the influence of Be star disk, it can not be explained 
within the inverse Compton scattering model of TeV emission from 
the system. 

Our XMM-Newton observations which were done during the 
period of first pulsar entrance to the disk (some 40 to 20 days 
before the periastron) revealed quite unexpected evolution of the 
spectrum of X-ray emission. In particular, the X-ray spectral slope 
hardened from r - 1.8 down to T ^ 1.2 roughly at the onset of 
the rapid growth of the X-ray luminosity. This behaviour can not 
be explained within the synchrotron model of X-ray emission from 
the system. We have shown that the observed hardening can be the 
result of the hardening of electron spectrum due to the Coulomb 
losses in the disk in the model where X-ray emission is the inverse 
Compton emission from the population of electrons responsible for 
the radio synchrotron emission. Otherwise, the hardening can be 
caused by the increase of the bremsstrahlung emission from these 
electrons. However, the bremsstrahlung problem faces a serious en- 
ergy balance problem. 
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